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Unconventional superconductors are most frequently identified by the observation of power-law 
behaviour on low-temperature thermodynamic or transport properties, such as specific heat. Here we 
show that, in addition to the usual point and line nodes, a much wider class of different nodal types 
can occur. These new types of nodes typically occur when there are transitions between different 
types of gap node topology, for example when point or line nodes first appear as a function of some 
physical parameter. We identify anomalous, non-integer thermodynamic power laws associated 
with these new nodal types, and give physical examples of superconductors in which they might be 
observed experimentally 



A defining feature of many unconventional Fermi su- 
perfluids and superconductors is the existence of lines or 
points on the Fermi surface where the energy gap van- 
ishes and so-called "nodal quasi-particles" can exist at 
arbitrarily low energies. In Fermi systems with such 
nodal quasi-particles the low temperature specific heat 
will show particular power law behaviours as a function 
of temperature. The expected point and line node power 
law dependencies were first derived in relation to the pro- 
posed low temperature superfluid phases of liquid '^He [T]. 
These were subsequently clarified further [5] and are now 
widely used to identify pairing states in unconventional 
superconductors Here we show that gap nodes in 

superconductors can occur in a number of more general 
types than simply the usual line or point zeros and that 
each of these has a corresponding thermodynamic sig- 
nature, typically in the form of non-integer power laws 
in low temperature specific heat. We predict that these 
anomalous power laws generically occur at points in the 
phase diagram where there is a topological change in the 
line or point nodal structure on the Fermi surface, which 
we illustrate with a specific example: the noncentrosym- 
metric superconductors Li2Pd3_:cPt2;B. In this case the 
gap node topological changes are also associated with 
changes of bulk topological quantum numbers for the 
quasi-particles. The future experimental observation of 
such non-integer power laws could therefore be used to 
identify not only superconductors with highly unconven- 
tional pairing symmetries, but also topological supercon- 
ducting and superfluid systems [6J. Similar physics are 
realised in the high-temperature cobalt doped Pnictide 
materials Ba(Fei_^T^)2As2 (T=Co,Ni,Pd) ^ZiiSj. 

For the usual types of point and line nodes the 
quasi-particle energy spectrum is linear in the vicinity 
of the gap node. A familiar example is the case of 
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the "^He A-phase with triplet pairing order parameter 
dk oc (fca; -|- iky)i, which has point nodes at the points 
k — (0,0, ±fci?) on the Fermi sphere, of radius kp- Near 
to these points the quasi-particle energy spectrum, i?k, 
obeys i?^ = Vp{kz — kpY + |dkP, where vp is the Fermi 
velocity, giving a linear dispersion relation as shown in 
Fig. [l]a). For this type of gap node we expect a specific 
heat, Ci, cx at temperatures T much lower than the 
critical temperature T^. The case of line nodes in the 
superconducting gap is also well known, as for example 
found in d-wave superconductors such as YBa2Cu307. 
Assuming in this case that the Fermi surface is an open 
cylinder of radius k p , and the energy gap is of the form 
Ak = Ao(A:^ — fc^), then there are four line nodes spaced 
around the Fermi surface at k = ^^-(±1, ±l)/-\/2. Again 
the energy spectrum is linear near these gap nodes, as 
shown in Fig. [l]c), and in this case we expect C„ oc 

While the generic point and line nodes of this type are 
expected in the majority of unconventional superconduc- 
tors, other types of nodal energy spectrum are possible. 
For example it is possible that the system has a point 
node, but that the energy spectrum is quadratic not lin- 
ear near to the nodal point, as shown in Fig. [l]b). The 
heat capacity at this critical point is a power law, with 
C„ (X rather than the usual oc for point nodes 
[5]. This aspect of non-linear excitation spectra was iden- 
tified by Leggett in his review of the superfiuid phases of 
■^He 12], though at the time no examples were known. 
Later a ground state with E2ti symmetry having this fea- 
ture was proposed for the heavy-fermion superconductor 
UPta [10] and possible thermal transport signatures were 
discussed [TP. Here we show that this different type of 
point node behaviour is not restricted to certain special 
symmetries but occurs generically when the gap struc- 
ture is changing from having a point node to having no 
nodes, as a function of some parameter defining the order 
parameter on the Fermi surface. We extend this notion 
to a much wider class of nodal types and node topology 
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FIG. 1: The four generic types of point and line nodes in a superconductor or Fermi superfluid and the two generic types 
of line crossings: (a) a point node with a linear excitation spectrum, (b) a shallow point node, (c) a line node with linear 
spectrum, (d) a shallow line node, (e) a crossing of ordinary hne nodes, and (f) a crossing of shallow line nodes. For the point 
and line nodes (a-d) we give example gap functions, densities of states and low-T specific heat capacity power law exponents n 
for Cv oc T". The coordinates kx\\ and ky\\ refer to a local coordinate system on the Fermi surface centered on the node, where 
kz is normal to the Fermi surface, and where k^^^ and fcyy are coordinates tangential to the plane of the Fermi surface at the 
node. The crossings (e,f) each have a distinct gap spectrum, density of states and power law which is distinct from the cases 
of point and line nodes alone, and which becomes dominant at low temperatures. The densities of states in these two cases are 
logarithmic, where the parameter L denotes the length of the line nodes on the Fermi surface (measured to midway between 
crossing points). These densities of states approximate a power law, as shown, and we also show the approximate heat capacity 
power law exponents n. 
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transitions. 

The fact that it is possible to continuously remove a 
point node of the gap is surprising, since the usual linear 
spectrum point node behaviour is usually thought to be 
topologically stable against perturbations \12^ . The shal- 
low node behaviour is not topologically protected, since 
it only occurs for specific parameter values in which ac- 
cidental cancellations occur between different order pa- 
rameter components. 

If we generalize the shallow point node behaviour to 
line nodes we obtain a spectrum near to the node shown 
m Fig. [l]d), which we will refer to as a shallow line node. 
In this case we find an anomalous non-integer power law 
for the low temperature specific heat, of the form C„ 
(X T^-^. Similarly to the shallow point node, this type 
of shallow line node only occurs at a boundary between 
a nodeless superconductor and one with ordinary line 
nodes, as discussed recenly for the pnictide supercon- 
ductors 13 IH]. In particular, both nodal and node- 
less compounds of the system Ba(Fei_a:Coj;)2As2 have 
been reported [T3l fT4]. suggesting that the T^^^ power 
law could be observed experimentally. If confirmed this 
would provide compelling evidence for the gap model 

lisj. 

Another type of nodal behaviour occurs when there 
is a crossing or topological reconnection of line nodes. 
These again lead to distinctive low temperature power 
laws, with non-integer power law exponents. A crossing 
of linear line nodes is illustrated in Fig. [l] (e). This type 
of crossing can occur in a d-wave superconductor, for ex- 
ample Ak oc k'^ — ky, in the case where the Fermi surface 
is spherical rather than cylindrical. Although the node 
crossing point occurs for a small portion of the Fermi sur- 
face, it turns out that this point dominates the low tem- 
perature heat capacity, and such crossings can therefore 
be detected experimentally. This is because the density of 
states resulting from the node crossing has a logarithmic 
term. This logarithmic term approximates a power law 
at low temperatures, and leads to a specific heat capacity 
contribution approximately of the form oc T^'*, which 
at low enough temperatures dominates over the usual line 
node contribution oc arising from the line nodes 
away from the crossing. Unlike the shallow point node, 
the crossing of line nodes is expected to be topologically 
protected against small perturbations. 

The final distinct type of gap node which we shall dis- 
cuss is a crossing of shallow line nodes, as shown in Fig. 
[T] (f ) . Again this leads to a logarithmic term in the den- 
sity of states, which dominates the low temperature be- 
haviour compared to the shallow line node case in the 
absence of crossings. The corresponding specific heat ca- 
pacity is approximately oc T^ '^ for the shallow line 
node crossing, compared to Cy oc T^'^ for the case with- 
out crossings. Interestingly, n = 1.4 is not so far from 
the value obtained for the ungapped Fermi surface of a 
normal (non-superconducting) Fermi liquid, n = 1. 

Fig. |2] shows the predicted power law temperature de- 
pendences in specific heat for each of the different nodal 
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FIG. 2: The simulated exponent of temperature T in the 
specific heat capacity as a function of temperature for each of 
the nodal characters. These match analytics (dotted lines) in 
the limit T — >■ 0. The methodology employed to produce this 
plot is described in detail in the SOM. 



types shown in Fig. [T] In all cases the predicted power 
law behaviour Cy oc T" is accurate at very low temper- 
atures, below O.lTc, and in many cases the low temper- 
ature power laws apply over a range of temperatures up 
to O.BTc. 

Non-integer power laws, especially powers below 2, are 
a highly distinctive signature of the low energy excitation 
spectrum, and their experimental observation would be 
a relatively clear and direct form of evidence for the ex- 
istence of the corresponding type of gap node in a given 
material. In systems where the gap topology is changing 
as a function of some experimental parameter (e.g. dop- 
ing or spin-orbit coupling) these anomalous powers will 
only be observed to the lowest temperatures if the ma- 
terial is exactly at the critical parameter value for the 
shallow node to occur. However, if the material is merely 
near to, not exactly at, the critical parameter value, then 
the shallow node power laws will hold over a wide temper- 
ature range, reverting to ordinary line, point, or nodeless 
behaviour at the very lowest temperatures. Again exper- 
imental observation of such a cross-over between power 
laws would be a clear signature of the unusual gap nodal 
structure - we provide an explicit example of this below. 

In addition to the ad hoc mechanisms mentioned above, 
we expect shallow line or point nodes to occur generically 
in superconductors with multi-component order param- 
eters. Consider the simple model shown in Fig. |3] This 
shows a hypothetical system, such as an ordinary s-wave 
BCS superconductor with a significant level of anisotropy 
in the gap. Here Aq is the average gap and Ai is a mea- 
sure of the gap variation with fcy parallel to the Fermi 
surface. In the case where |Ao| > |Ai| the gap is node- 
less everywhere on the Fermi surface. On the other hand. 
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FIG. 3: The transition between a nodal and non-nodal gap 
structure in any superconductor as a result of continuously 
admixed order parameters will be characterized by a non- 
linear gap energy spectrum close to the node. This is because 
in any wave-like order parameter the deep minima have no 
linear order term in the expansion. 

if |Ao| < I All then the gap changes sign at some points 
on the Fermi surface, leading to pairs of line nodes. At 
the critical parameter value |Ao| = |Ai| the line node 
pairs coalesce, leading to a single, shallow, line node. In 
this example no symmetry change occurs in the gap func- 
tion, however "accidental" line nodes occur within a single 
symmetry of the Cooper pair order parameter, because 
of the large degree of gap anisotropy. 

An explicit example of this mechanism by which 
anomalous power laws may be observed is provided by 
non-centrosymmetric superconductors where the super- 
conductor order parameter 

Ak = (Ao,k + dk • cr) lo-y (1) 

has singlet and triplet components Ao,k and dk, respec- 
tively, such as the compounds Li2Pd3_a;Pta;B |16fll8| . 
Platinum doping in this system is thought to increase 
the spin-orbit coupling, and hence increase the triplet 
component of the order parameter relative to the singlet 
component |ilEM21j . In doing so the spin-orbit cou- 
pling introduces nodes in the gap. Within a standard 
symmetry classification of the possible order parameter 
structures the simplest order parameter takes the form 

[23121] 

Ao,k = Ao (2) 
dk = Ao{A[X,r,Z]- (3) 
B [X {Y^ + Z^) ,Y {x^ + z^) ,z (y2 + x'^)] } 

where X, Y, Z are functions with the same symmetry 
properties as the components of the electron wave vector 
kx,ky,kz and A and B are material-specific parameters 
that depend on doping x. A large number of different 
gap nodal structures are obtained, as shown in Fig.|4](a). 
As a function of the parameters A, B we obtain the rich 
phase diagram of Fig. |4] (b). The corresponding specific 
heat exponents are illustrated in Fig. |4] (c) . Upon dop- 
ing, the material explores the phase diagram along some 
unspecified path from small values of A and B (x — 0) to 
larger values (x — 3). If the triplet gap components are 
small we have a nodeless but anisotropic superconductor. 

Inspection of the phase diagram shows that at some 
values of A and B we expect point nodes to appear, la- 
belled by T in Fig. |4] (a,b), before line nodes appear 



on the gap, labelled by '2' in the same figure. These 
point nodes are shallow point nodes, and so instead of 
the usual power law dependence of the heat capac- 
ity for point nodes, we instead expect a power law 
dependence. Upon further change of the A and B pa- 
rameters these shallow point nodes expand into rings of 
ordinary line nodes, as shown in Fig. [4](b), plot labelled 
'2'. Increasing the parameters yet again there is a level 
which causes a topological line reconnection, labelled '3', 
where the different line nodes intersect before becoming 
rings again, now centred at different points on the Fermi 
surface, '4'. Note that the number of line nodes on the 
Fermi surface is always a multiple of two as predicted by 
the topological theory of gapless phases in time-reversal- 
invariant superconductors [22'. That number jumps at 
the line reconnection transition, indicating its topologi- 
cally non-trivial nature. The line reconnection transition 
has a distinct thermal signature in the low T heat capac- 
ity originating from the crossing of line nodes associated 
with such a line reconnection. As is evident from the 
figure, there are many other paths along the phase dia- 
gram with similar phase transitions. In our analysis of 
the cubic non-centrosymmetric superconductors we also 
found a single instance of crossing of shallow line nodes, 
labelled '11' in Fig.gl These are realised when line nodes 
are spontaneously introduced on a nodeless section Fermi 
surface without first appearing as points. The specific 
heat exponent of n = 1.4 for this case is clearly seen in 
Fig.|4](c). 

Schnyder et al. |23l have shown that a number 
of these phases of cubic noncentrosymmetric supercon- 
ductors possess non-trivial topological quantum numbers 
[HI [22] 1 and have corresponding surface edge states of 
topological nature. In principle, these surface states can 
be probed experimentally by electron tunnelling through 
an interface with a normal metal [25j. However, we pro- 
pose that a more direct, bulk, measurement of this non- 
trivial gap topology can be obtained by measurements 
of the low temperature heat capacity. If the anomalous 
power laws we have predicted can be observed, then it 
should be possible to detect the phase boundaries be- 
tween different types of topological gap structures, and 
hence determine in a bulk measurement the topologi- 
cal quantum numbers associated with different doping 
regimes. Phases 'l'-'5' in F ig. [4] co rrespond to those con- 
sidered by Schnyder et al. [23|l24]. 

The anomalous power-laws due to changes of nodal 
topology are realised exactly at each transition. On the 
other hand, by being close to, but not exactly at one 
of those transitions one expects a crossover from con- 
ventional (fully gapped or normal power-law) behaviour 
to anomalous behaviour as the temperature is increased. 
A demonstration of this is provided for the case of 
Li2Pd3„2;Pt2:B -see Fig.|4](c), specifically the curves cor- 
responding to A — 3, B = 3.5 — 4.0. This may be used 
to detect topological transitions in superconductors. The 
detection of the phase boundary itself is facilitated by its 
higher specific heat, at fixed low temperature, due to the 
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FIG. 4: (a) Fifteen gap topologies corresponding to different triplet admixtures in a cubic non-centrsymmetric superconductor. 
Within this set the only point nodes are shallow point nodes, occuring in plots 1,5 10 and 15. Topologies 3,6,8 and 14 are all 
line reconnection transitions between distinct topologies and include crossings of linear line nodes. Topology 11 is a crossing of 
shallow line nodes, (b) Phase diagram of gap node topology corresponding to triplet admixture parameters A and B in Eq. (|3| 
on a spherical Fermi surface. Each line represents a phase boundary where gap nodes appear or disappear on the Fermi surface, 
(c) The specific heat exponent as a function of temperature for different values of the A and B parameters, as indicated. The 
inset shows the value of the specific heat as a function of B for fixed A = 3 and constant temperature, T — 0.005 Tc- Details 
of the calculation of the quasiparticle spectrum and phase diagram are given in the SOM. 
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enhanced phase space available to quasi-particles in that 
state. This is illustrated in the inset to Fig. |4] (c). We 
emphasize, however, that it is in the anomalous power- 
laws that we find the evidence of the special nature of 
quasi-particles at the phase boundary. 

In recent years it has become widely recognised that 
point or line nodes in the bulk spectrum of topological 
matter may evolve into topolgically-protected surface or 
vortex core states with zero dispersion |26l l?7| . In con- 
trast, the shallow dispersions discussed here exist in the 
bulk excitation spectrum near topological transitions. 

To conclude, we have shown that a more general clas- 
sification of nodal characters of superconductors yields 
anomalous, non-integer power laws for low-temperature 
thermodynamic quantities. Some such power laws can be 
used to detect line node crossings; others are expected at 
phase transitions where the nodal topology reconfigures 
itself. An example is provided by the topological transi- 
tions expected to occur as a function of doping in the non- 
centrosymmetric cubic superconductor Li2Pd3_a;Pta;B. It 
is important to stress that the observation of anomalous 
low-temperature power laws, just as the ordinary, integer 
ones, requires going to temperatures well below Tc- On 
the other hand, it is not necessary to sit exactly at the 
transition, as there is a crossover temperature separating 
conventional from anomalous power-law behaviour that 
converges to zero there. In that sense the influence of the 
topological transition extends throughout the bulk phase 
diagram in a manner similar to that of a quantum critical 
endpoint [28,,29j. 
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Supplementary Material 

I. CALCULATION OF SPECIFIC HEAT OF 
NODAL SUPERCONDUCTORS 

To describe the energy Ei^ of quasi-particles near any 
specific point on the Fermi surface it is convenient to 




where /3 — l/ksT. The low T specific heat capacity can 



change coordinate systems to three mutually orthogonal 
k-vector components: one perpendicular to the Fermi 
surface, and two parallel to it fej-y, We then have 

-^k — vp'^k^^ + IS.\^ where vp is the Fermi velocity, fcj^ — 
when |k| — kp (the Fermi vector), and is a function 
determined by the superconducting gap function, which 
may be either singlet, triplet or a combination (for non- 
centrosummetric systems) . 



The gap energy spectra shown in Fig. 1 (main text) are 
described by an equation for the gap energy Ak which has 
the form A^ = fcj^n + ... near to a gap nodal point, 

where /^i, is a positive tensor of rank 2 which is deter- 
mined by the symmetry of the gap function around that 
nodal point. In general /^i, will have two principal axes, 
which we can choose as fc^y and fcyy. If the two eigen- 
values corresponding to these axes, Ix and ly are both 
positive and non-zero we have an ordinary point node 
with linear spectrum illustrated in Fig 1(a) (main text) 
(with /i = \/IxIy) ■ If one eigenvalue is positive and non- 
zero (/i) and the other is zero we have an ordinary line 
node. In the case of two zero eigenvalues it is necessary 
to continue the exansion to quartic or higher order, giv- 
ing the various shallow nodes or crossings illustrated in 
Fig. l(main text) . In principle even higher powers may 
be possible, but we have not found any realistic examples 
of gap models where this could occur. In contrast, all the 
node types presented in Fig. 1 (main text) can be realised 
in the specific gap models presented in this paper. 



Computation of the density of states integral g{E) = 
/ / j 6{E\^ — E)dkxdkydkz for each of these gap energy 
spectra gives a unique expression. The specific heat ca- 
pacity in the superconducting state is in turn given by 

El 



i?ksech^(^) (4) 
I 

therefore be found analytically from g{E) by evaluating 
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the integral ifcs/J^ dE g{E)E'^ sech^ (^£^/2) , giving 
the power law dependence of a low T specific heat ca- 
pacity measurement. We find the terms listed in Fig. 1 
(main text) for point and line nodes with linear and non- 
linear excitation spectra, where L is the length of the line 
node from end to end. For cases where line nodes cross we 
expect the terms listed in Fig. l(main text), where L is 
now the length of the line node on the Fermi surface mea- 
sured from the centre to the corner of the crossing whose 
contribution is being calculated. The cases in Figs. 1 
e) and f) (main text) are not power laws, however they 
approximate very closely to the non-integer power laws 
indicated, with very slight variations depending on the 
constant L. In the case of line reconnection transitions 
such as Fig. 4, plots 3, 8, and 14 (main text), where 
the Fermi surface includes regions where the gap energy 
spectrum expands as — Iikx\\^ky\\^ , the contribution 
of the line crossing to the density of states is ~ E^-^ . 

The low-T specific heat capacity can also be calcu- 
lated numerically using Eq. ^ by summing over a suffi- 
cient number of integration points corresponding to dif- 
ferent values of the gap in k-space near the Fermi surface 
for selected temperature ranges below Tc. The expo- 
nent of temperature, n, can then be obtained using the 
general formula n = dlnCu/dlnT. When this quantity 
approaches a constant we have cx T". This yields 
the curves in Figs. 2 and 4 (c) (main text), where a 

1/2 

Ginzburg-Landau temperature dependence, (T — Tc) , 
was assumed for the overall amplitude of the gap energy. 
This introduces only a weak temperature-dependence for 
temperatures T <^ Tc- 

Each contour in Fig. 4 (a) (main text) corresponds to a 
line where the gap energy Ak derived from the gap func- 
tion of Eqs. (2,3) (main text) vanishes. This gap energy 
is given by |Ak| = |Ao ± d(k)| = 0, where d(k) = |dk|. 
The phase diagram in Fig. 4 (b) (main text) shows the 
values of the parameters A and B where the nodal topol- 
ogy changes for this form of the gap energy. This form 
of Ak has also been used to compute the specific heat 
exponent in Fig. 4 (c) (main text). 



II. SPECTRUM OF O-GROUP 
SUPERCONDUCTORS 

Here we give the details of the derivation of the nodal 
topologies and phase diagram presented in Fig. 4 (b) of 
the main text. We use a single band mean-field model 
with antisymmetric spin-orbit coupling to obtain the re- 
sulting gap node topologies on the Fermi surface, which 
we take to be spherical for simplicity. The Bogoliubov-de 
Gennes Hamiltionian is 



- ( Mk) A(k) \ 

^("^^ - I^At(k) -r(-k)j- 



The noncentrosymmetric (NCS) gap function which 
mixes singlet and triplet as a result of the parity violation 
in the crystal lattice structure is given in Eq. (1) of the 
main text, where <yx,y,z denotes the Pauli matrices and 
Ao,k and dk denote the singlet and triplet components 
of the gap respectively. 

The normal state Hamiltonian h(k) given in [53' for 
non-interacting electrons in a crystal without inversion 
center is 

h(k) = £k/ + 7k • (6) 

where ek = £-k is the metallic dispersion energy. The 
second term represents the antisymmetric SOC where the 
coupling constant 7k = ~7-k- We will assume the d- 
vector dk to be parallel to 7k. This ensures that SOC is 
not destructive to the triplet component of the NCS gap 
function f23l|30]. 

In a general Ginzburg-Landau Theory the possible 
symmetries of the superconducting instability depend on 
the symmetry of the normal state. In this case the sym- 
metry of the normal state is given hy Gn = OxTx U{1) 
where O denotes the NCS cubic space group, T the time 
reversal symmetry and U{1) the gauge symmetry oper- 
ation. In any superconducting transition the U{1) sym- 
metry is broken however other symmetries of the normal 
state may be broken as well. For the purposes of our anal- 
ysis we will, following look at the case where only 
the gauge symmetry and no further symmetries of the 
crystal lattice are broken i.e. the superconducting insta- 
bility has the symmetry of the Ai irrep of the O group. 
This means that the gap function, given in general by 
equation (1) of the main text, takes the specific form of 
Eqs. (2,3) of the main text. 

With dk and 7k explicitly defined the Bogoliubov de- 
Gennes (BdG) Hamiltonian can be diagonalized. The 
single band BdG hamiltonian with SOC defined parallel 
to the d(k) component of the gap ^23. .30] has the prop- 
erty that splitting of the Fermi surface occurs as a result 
of SOC whereas zero energy excitations occur as a re- 
sult of the combination of singlet and triplet components 
of the gap function. In this case the diagonalized BdG 
Hamiltonian has eigenvalues of the form 

i^ = ±\/(£k-A^±l7kl)' + |Ao±|dkir (7) 

The SOC term 7^. in this expression has the effect of split- 
ting Fermi sheets whereas the second term |Ao± \dk\ \ de- 
termines the node topology of the superconducting gap. 

The zeros of the term |Ao ± |(ik|| are the locii of zero 
energy excitations on the Fermi surface. In order to en- 
sure that the SOC parameter does not significantly warp 
the spherical Fermi surface and therefore change the node 
topology we assume that ek ^ 7k ^ Aq. This hierar- 
chy of energy scales is reasonable to assume as ek is of 
the order of eV, Aq is of the order of meV, and observed 
SOC-induced Fermi surface splittings are small relativis- 
tic corrections to ek- 
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The conditions for zeros in the gap energy spectrum 
are Aq — |dk| = and Aq + |dk| = 0. These conditions 
produce equivalent phase boundaries in the phase dia- 
gram. The gap function is then normalized about the 
singlet gap: 



text separate regions of constant gap topology which were 
obtained for critical values of the polar and azimuthal an- 
gles corresponding to the points where nodes first appear 
on the gap or where line nodes reconnect at a topological 
transition. 



^0 



1 - 



Ao 



= A2|i-|d;,ll^ = o 



(8) 



The term has here been replaced with the term |dj^| 
whose coefficients are scaled by a factor of Aq. This 
expression can be fully expanded and used to derive 
the topological phase boundaries by substituting in the 
vector d'(fc) = [AX - BX{Y^ + Z'^),AY - BY{Z^ + 
X'^),AZ - BZ{X'^ + y^)]- A coordinate transformation 
to spherical polar coordinates using X — cos((/)) sin(6'), 
Y = sin((/)) sin(6') and Z — cos{9) allows us to select spe- 
cific points on the Fermi surface and find solutions in A 
and B. The phase boundaries of Fig. 4 (b) of the main 



Phases 1-5 in Fig. 4 (a) of the main text correspond 
to those considered by Schnyder et al. |23l [21] . The 
coordinates 9 = cos^^ (I/a/S) ,4'=j correspond to one 
topologically distinct set of point nodes on the spherical 
Fermi surface shown in Fig. 4 (a), plot 5 of the main 
article. Substituting these into the equation Ao — |dk| =0 
gives the relationship between coefficients 'A' and 'B': 
A — ± 1. The line reconnection shown in Fig. 4 (a), 
plot 3 of the main article corresponds to the coordinates 
9 — I , (/) = and produces solutions A — ^^^^ . The final 
phase boundary in Fig. 4 (b) corresponds to point nodes 
at the coordinates 9 = |, (/> = shown in Fig. 4 (a), plot 
1, which produces the solutions A = ±1. 
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